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In this paper, we study 3-point correlation function of primordial gravitational waves generated 
in the de Sitter background in the framework of the general covariant Hoi^ava-Lifshitz gravity with 
an arbitrary coupling constant A. We find that, to the leading order, the interaction Hamiltonian 
receives contributions from four terms built of the 3-dimensional Ricci tensor Flij of the leaves t = 
constant. In particular, the 3D Ricci scalar 7? yields the same fc-dependence as that in general 
relativity, but with different magnitude due to coupling with the U{1) field A and a UV history. 
The two terms RijR^^ and (Vi?^*") (ViRjk) exhibit a peak at the squeezed hmit. The term R)R?^R'l 
favors the equilateral shape when spins of the three tensor fields are the same but peaks in between 
the equilateral and squeezed limits when spins are mixed. This is due to the effects of the polarization 
tensors. Hence, a detection of the squeezed shape or absence of in the bispectrum cannot rule out 
completely higher derivative gravity theories, at least for R^ theories. The consistency with the 
recently-released Planck observations on non-Gaussianity is also discussed. 



I. INTRODUCTION 

The Hofava-Lifshitz (HL) theory of quantum gravity, 
proposed recently by Hofava 1 , motivated by the Lif- 
shitz scalar field theory in solid state physics [5], has 
attracted a great deal of attention, due to its several re- 
markable features in cosmology as well as some challeng- 
ing questions, such as instability, ghost, strong coupling, 
and different speeds in the gravitational sector [5]. To 
resolve these issues, various models have been proposed, 
along two different lines, one with the projectability con- 
dition gHH], 



related to them, including instability, ghost, and strong 
coupling, are resolved automatically. This was initially 
done with A = 1, where A characterizes the deviation of 
the theory from general relativity (GR) in the infrared, 
as one can see from Eqs.(2.1) and (2.2 1 given below. It 



N = N{t), 



(1.1) 



and the other without it [5l412j. where N denotes the 
lapse function in the Arnowitt, Deser and Misner de- 
compositions [13 . In particular, Hofava and Melby- 
Thompson (HMT) proposed to enlarge the foliation- 
preserving diffeomorphisms of the original model. 



6t = -f{t), 6x' = -Cit,^), 



(1.2) 



was soon generalized to the case with any A [7| , in which 
it was shown that the spin-0 gravitons are also elimi- 
nated fTlIH], so that the above mentioned problems are re- 
solved in the gravitational sector even with any A. In the 
matter sector, the strong coupling problem, first noted 
in [5], can be solved by introducing a mass Af* so that 
M^. < A^, where M^, denotes the suppression energy of 
high order operators, and A^ the would-be energy scale, 
above which matter becomes strongly coupled [Tlj, sim- 
ilar to the non-projectability case without the enlarged 
symmetry [T5]. The consistence of this model with so- 
lar system tests was investigated recently [16] , and found 
that it is consistent with observations, provided that the 
gauge field and Newtonian prepotential are part of the 
metric, so that the line element ds^ is invariant not only 
under the coordinate transformations (1.2), but also un- 



often denoted by Diff(M, J^), to include a local U(l) 
symmetry, so that the reformulated theory has the sym- 
metry [5], 



f/(l) K Diff(M, J"). 



(1.3) 



With such an enlarged symmetry, the spin-0 gravitons, 
which appear in the original model of the HL theory [1] , 
are eliminated [5j [6], and as a result, all the problems 
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der the local U(l) gauge transformations, 

S^N = 0, Sa^N, = NV,a, S^g,^ = 0, 

5aA = a - TV^V^a, Saip = -a, (1.4) 

where a denotes the U(l) generator, and d = da/dt, 
and N'^ and gij are, respectively, the shift vector and the 
3-metric of the leaves t = constant, with Ni = gijN^ . 
A denotes the U(l) gauge field, and Lp the Newtonian 
pre-potential. For detail, see [H]. 

In this paper, we shall work within the HMT frame- 
work of the HL theory with the projectability condi- 
tion [5HH], although our main results are expected to 
hold equally in the non-projectable case [TU], as the 
tensor perturbations are almost the same in both cases 
O [ini HI] . For the current status of the models proposed 
in H [i m H Uni , we refer readers to [M [HI ■ 



Continuing our previous study of the statistics of the 
primordial perturbations in single-field slow-roll inflation 
in the HMT framework [TSl[Tn], we study here the 3-point 
function of the tensor perturbations in de Sitter back- 
ground. Non-Gaussianity of tensor perturbations have 
been studied intensively recently in various theories of 
gravity. But, as far as we know, this is the first time 
to investigate this problem within the framework of the 
HL theory. In particular, [20] studied the non-Gaussian 
characteristics in the most general single-field inflation 
model with second-order field equations and found that 
the interaction at cubic order is composed of only two 
terms, which generate squeezed and equilateral shapes. 
On the other hand, [21] focused on possible parity viola- 
tions. In this paper we shall focus on the shapes of the 
bispectrum generated by the various terms in the gravity 
sector, especially the higher spatial derivative terms and 
leave the topic of parity violations to future studies. For 
detail on the power spectrum and its scale dependence of 
the tensor perturbations, we refer readers to section 6.2 
of [H]. 

The rest of the paper is organized as follows. In Sec- 
tion II we give a brief review of the non-relativistic gen- 
eral covariant Hofava-Lifshitz theory of gravity with the 
projectability condition and an arbitrary coupling con- 
stant A. The interaction Hamiltonian Hi is analyzed in 
Section III, and was found to receive contributions from 
only four terms in the potential part of the theory. Sec- 
tion IV performs the integration of the mode functions 
of the quantized fields, where we also present conditions 
with which the UV history in the integral can be ignored 
as a subleasing error term. We then plot various shapes 
of the bispectrum generated by the four terms in Sec- 
tion V. In Section VI, we consider the non-Gaussianity 
constraints from Planck observations released recently to 
obtain restrictions on the energy scale M* and the infla- 
tion energy H. Finally, in Section VII we summarize our 
main results. 



II. GENERAL COVARIANT HL GRAVITY 
WITH PROJECTABILITY CONDITION 



The action of the general covariant HL theory of grav- 
ity with the projectability condition can be written as 

[SHa, 



S = -^^ I dtd'xN^^CK-Cv+C^ + CA + C 

dtd^xN^CM: (2.1) 



where g = det{gij), G is the Newtonian constant, and 



Cx 



C'.9o + giR + ^ (52 i?' + 93R^JR'') 

+4 (ff4i?' + 55i? R^JR'' + geR^RiRi 
+ ^ [grRAR + gs {\7,R,k) {S/'W")] , 



Ca - ^{2^9-R 

Cx = {l-X)\{Aipf + 2KAip 



(2.2) 



(2.3) 



Here A characterizes the deviation of the theory from 
general relativity (GR) in the infrared, as mentioned pre- 
viously, A = g'^ViVj, Ag is a couphng constant, the 
Ricci and Riemann tensors Rij and i?*^./ all refer to the 
3-metric gij, and 



K,, 



2N 



- -gijR + Aggij. 



The coupling constants gs (s — 0,1,2, 
mensionless, and 



A 



90 



(2.4) 
>) are all di- 

32.G' (2-^) 

is the cosmological constant. The relativistic limit in the 
IR, on the other hand, requires, 

31 = -1. (2.6) 

Cm is the Lagrangian of matter fields, and for a scalar 
field X, it is given bv [T81I22]. 



Cm = 4°) +4^-^), 



£?) 



2iV2 



(x-iv*va) -V, 



V = Vix)+i- + Vi (x) j (Vx)' + V2 (x) vf 
+V3 (x) T^f + Vi (x) V2 + V5 (x) (Vx)'7'2 



-V6{x)ViV2, 



(2.7) 






A- A 

N 
_l 

N 

4 



ci(x)Ax + C2(x)(Vx)' 



(x"A^'V,x)(vV)(Vfcx) 



^V'^»(Vfex) 



(2., 



with ci (x) , C2 (x) , y{x) E^nd Vn{x) being arbitrary func- 
tions of Xi and 

Vn = A"x- (2.9) 

For detail, we refer readers to [H [18] . In the following 
we shall use directly the symbols and conversions from 
these papers, without further explanations. 



III. THE INTERACTION HAMILTONIAN 

We assume that no tensor perturbations exist in the 
matter sector, and the tensor perturbations of the metric 
around a spatially flat FLRW metric is [20] 



N 



N' 



0, gij 






where 



h,.-i + -h^kh'^ + -hijW^hki + 



i-ij 



(3.1) 



(3.2) 



The metric perturbation is defined in this way such that 
there is no cubic term involving two time derivatives in 
GR [53]. The small dimensionless quantity hij satisfies 
the transverse and traceless condition d^hij — = h]. 
Moreover, its indices are lowered and raised by dij and 
d''^ . Thus for simplicity of notations, in this paper we 
shall not distinguish super-indices and sub-indices, with 
the understanding that when an index appears twice, 
summation over that index is performed. We further 
introduce a short hand notation 



(^"')ii = ^ifci^feifes • • • ^fc„-ij 



(3.3) 



With these perturbations, the interaction Hamiltonian 
at leading order is found to receive contributions from 
only four terms 

R, R,jW\ R]RiRt {V'W''){V^RJk)■ (3.4) 

The kinetic part of the action Sk, like the case in GR, 
does not contribute to Hi even though A now could differ 
from 1 [cf. Eq.(IZ2])]. We find that. 



R 



-h'^'''''{2h,ahjb-h,,hab).. 



(3.5) 



R^jR'^ = a-4 



{d^h'^] 



, '^ia.h'^jh.a 



'ab 



^ia,jb 



.iT'ab,ij 



rt ^i'j.ab 

2 



(3.6) 



R'jRiR., 



[dHh,,)][d'{h,k)][dHhk.)],{3.7) 



(V'i?^"'=)(V,i?,fc) = — (d'h,,) [h'^" (d'h^j)^^, 



Ah' 



'■" {d^h 



bj) 



{d% 



^ia^h'^jh.a 



in an essential way. The model considered in [^D] de- 
scribes the class of single-field inflation models where the 
coupling between the inflaton and gravity could be differ- 
ent from the canonical form while Lorentz symmetry is 
kept. This explains why their interaction term posses 
time-derivatives. On the other hand, coupling of the 
gravitational sector with scalar matter in our model is 
in the canonical form [cf. Eq. (2.7)]. The higher spatial 



derivative terms exist because of the Lorentz symmetry 
breaking. 

We define the Fourier image of hij here in the canonical 
form, 



'^-■-/(Sa E_ 



mhlit)e 



^kx 



(3.9) 



where h^(t) is now a scalar quantity, and the rank-2 
polarization tensors satisfy ef,i(k) — fc'ef (k) = and 



4(k)£f,(k) 



} By a proper choice of phase, they 



satisfy the relation ^20j 

[4(k)]* = £^^(k) = 4(-k). 

Then, in the momentum space we find 



(3.11) 



a^R-^ — 



a^W^R, 



Ci 



(123) - 
q 4 



Co 



(123) 



C. 



(123) ^-h cyclic, (3.12) 



9i + 92 + ll 



6 



Ci 



(123) 



1 r 
'4 



a RjRj^Ri 



C2 

-\- cyclic, 

(919293) 
24 
+ cychc. 



(123) 



C3 



(123) 



C. 



(123) 



(3.13) 



(3.14) 



{V'W') (V,i?,fe) ■ 



4 I 4 I 4 
9i + 92 + 93 



Ci 



(123) 



Co 



+ 



9 ? 

9i93 
6 

9 9 

9i92 



(123) - - 
q 4 



C. 



(123) 



6 

-|- cyclic 



(123) - - 
q 4 

(123) - I 
q 4 



C3 
C2 



(123) 



^(123)| 



(3.15) 



-h""' 



'^^ia,bj I ^'^ij^ab v i^ab^ij 



(3.8) 



The contribution from R is the same as that obtained 
1 [2^. The contribution from R^.R-j^R^ differs from theirs 



^ This set of polarization tensors is related to those introduced in 
[T8| through 

4 = 2(4^-5)- (3-10) 



where "cyclic" refers to the cychc rotation of (1, 2, 3), 
and we've defined the symbol -^^ to have meaning of 



3 

n 









(3.16) 



and introduced the shorthand notations 






(123) ^ fc,„fc,,£5(qO£S(q.)e;g(q3), 

^(123) EE A:,„fc,,£5(q,)<Uq2)<|(q3), 

(123) ^ fc,„fc,,£5(q0e5(q.)<Uq3), 

(123) ^ A;,.fc,,£5 (q,)<Uq,)<Uq3) , 

(123) ^ e5(qO£J|(q.)eS(q3). 



(3.17) 



Hence, the Hamiltonian reads 

(27r)3j(q, +q, + q3) 



Hi{t') = 



a-^C 



3/^2 



/nl:7^^«(.') 



.93 



(A-l)i?+p(ff^i?, 



+ ^(v^ff'=)(v.i?,,) + ^i?;.i?ii?: 



56 nj nj pfe 



C' 



(3.18) 

where expressions of R, [W^ R jj) , { \7^W''){ViRjk) and 

R}RiRf are given in Eqs.( |3.12ptl5| ). 

Promoting the scalar variable ft^(tj to a quantized field 



hUt) = /ik(t)a.(k) + r k(i)al(-k), 



(3.19) 



where the creation and annihilation operators satisfy the 
commutation relation 



a,(k),4(k')l ^{27rfS{k-k')Sss', 



(3.20) 



we are now in the position to employ the in-in formalism 
[23] . The 3-point correlator we seek can be calculated at 
leading order 



I I dt' 



i(27r)^(5(ki 



^^.)jy^'^ 



e 



W{K,t'-t)-W* } 



x{F,,,,,3(-K,f') 

+ 5 permutations of (ki,k2, ka), (3.21) 
where we have defined the product of mode functions 



W{K,t';t) EE hk,(t')hl^{t)hkS')K.(t)hkAt')hl^{t), 

(3.22) 



and for the contraction {(qi,ki), (q2,k2), (qsjka)} 



-^^3132^3 \ -I^I ' ) 



a?{t 



1\ 1 
6 



Cl - -C2 - 1^3 



(123) 



33 kl + kl + kl 
a4(t')C^ 6 

58 fcf + fc| + kt 



6 



2f,2 



9s kfk. 
' a^{t')C^ 6 

9% klkl 
' a^{t')C'^ 6 

ffg kl K2 k^ 

'a^{t')C^ 24 



Cl - -C2 - -C3 

.^1 " 4^2 - ^C3 



(123) 



(123) 



Cl - -C2 



Cl - -C3 



(123) 



(123) 



C. 



(123) + cyclic, 



(3.23) 



with 



Cl 



J123) EE fc,„fc,,£5(-k0£^^(-k,)ejg(-k3) 
= fc,„fc,,£,::.^Mk,)£r/Mk.)e-,^Mk3) 



— kiakit 



£5(k0<^(k,)e^^g(k3) ,(3.24) 



and similarly for 



C2 



C3 



(123). We see 



(123) and 

that the first line of ( |3.23| is contributed by R, the sec- 
ond line by R^^Rij, the third, fourth and fifth lines by 
(yiRJk-^ {S/iRjk), and the last line by R^RlRt 



IV. THE MODE INTEGRATION 

As was noted in [12], the /c-dependence of the bis- 
pectrum receives contributions from both the interaction 
Hamiltonian and the mode function, whose behavior is 
determined by the EoM and its initial state (mathemat- 
ically the initial condition). In principle, one should per- 
form the integration in (3.21 1 tracing the behavior of the 



mode function through out its history from the initial 
time (ti) until several e- folds after horizon exit when the 
mode freezes {t). However, during the early times the 
modes are highly oscillatory and the integration during 
that time should be mostly canceled out by itself. This, 
of course, requires some conditions on the integrands. 

Looking at the integrands in ( 3.21[ ), we see that the 
non-oscillating time-dependent parts come from the com- 
mon factor a^{t'), the factor of a~'^"{t') in Fs-^S2S3 (n — 
1, 2, 3) and from the three mode functions. The oscillat- 
ing time-dependent part in the mode function is modeled 

as [ig 



ii„i\i} 



h(x{r]')"'exp[iXik\r]') 



(4.1) 



where m and / are dynamical, i.e. they change when 
different terms dominate the dispersion relation. Since 



we are working with de Sitter space, then a = —l/{Hri), 
hence the integration can be written in a schematic form 



drj' exp[iXi [ki 



Ih l3rn+{2n-4) ^ 



(4.2) 



When 3m + (2n — 4) < 0, the integration must be small 
because yy'3™+(2"-4) changes very slowly during the early 
times {ri' -^ — cxo) while the oscillating part has high fre- 
quencies; similarly when < 3m + (2n — 4) < I, the oscil- 
lation smoothes out the changes in 7y'3™+(2n-4) g^g Q^g (^^-^ 
always perform a change of variables r = 7y'3m-i-(2n-4)+i 
and integrate over r with a purely oscillating integrand. 
When 3to -|- (2n — 4) > I, on the other hand, care is 
needed. However, if 3to -|- (2n — 4) is not so higher than 
I, the errors introduced by ignoring this history should 
not be too large. In fact for the current model, when 
k^ term dominates the dispersion relation, to = 0, ^ = 3 
and n„Ax = 3, hence 3m -I- {2n — 4) < 2 < L Therefore 
we shall work under this assumption below, and consider 
only the period when the fc^ term dominates the disper- 
sion. During this period, the oscillating mode function 
takes the form. 



Vkit) = 



C+{kv,) 



V2k 
C^ikju 
2k 



1- 
1 



1 



ic-j-kr] 
1 



icTkf] 



'icT^rj 



JCTk'n 



(4.3) 



where c^ = (1 — A), the canonically normalized field Vk{t) 
is related to hk{t) through^ 



Vkit) 



aM, 



pi 



hk{t), 



(4.4) 



and the constants C+ and C_ are in general functions of 
iJ, Mpi, the new mass scale M* and transition time rji 
[l9] . The reason for this dependence on {kiji) is the UV 
stage when the dispersion relation differs from the rel- 
ativistic form significantly. The detailed form however, 
depends on the procedure of matching of this "relativis- 
tic solution" with the solution in the UV region. One 
example is the matching we considered in [19] . 

Summing up the discussions, we see that the UV his- 
tory has at least three effects on the integration we are 
considering: one is that rji can no longer be extended to 
Euclidean space — oo(l -I- ie); second, the dependence of 
the normalization on rji] third, if the mode underwent 
a non-adiabatic period (o;^ < 2/77^), then a "negative 
frequency" branch exists (C_ 7^ 0). Since we've seen in 
[in] that a mixture of "negative frequency" and "positive 
frequency" modes gives an enhanced folded shape in the 
bispectrum |19J , in this paper we focus on the case when 
C_ =0. 



^ Note that the relation here is different from the relation in 1181 
due to the different normalization we choose for the polarization 
tensors. The EoM's however, are the same for both. 



The bispectrum is then given as 



{K\it)Klit)hg{t)) 

{2TTfS{k,+k,+k,)C 



-£(UV) + f (Finite rji) -t- 5 perm.'s. 



V2HC^ 



2G 



S1S2S3 



(K) 



ht<J htO httj 

"'1"'2"'3 



(4.5) 



where we have introduced two error terms signifying 
that we considered only the relativistic region in the 
integration, and defined, for the contraction pairing 
{(qi,ki),(q2,k2),(q3,k3)}. 






(K) 



(^-1)^ 



Ci - -C2 - -C3 



(123) 






,IZ 



,IZZ 



nJ-j ~| hjo 



Ci - -C2 - -C3 
Ci - -C2 - -C3 



(123) 



(123) 



C 6 

"'1"'2 



Ci - -C2 



\xxx 



i.2pjN2 



Ci - -C3 



(123) 



(123) 



21,2 7,2 r 



—,XX1- 

C* 24 



G. 



(123) -I- cychc. 



(4.6) 



where 



XX = -2 



XXX = 8 



- (fci -I- fc2 -t- fca) + fclfc2fc3 
Cr^ (fci + fc2 -I- fes) 

(fci + A;2 + k-i) {k\k2 -f fc2fc3 + ^3^1) 

Cr^ (fcl -I- ^2 + kz) 
(fci + fc2 + fcs) +3/i;ifc2fc3 

Cy (fci + fc2 + fcs) 
(fci -I- fc2 -t- fc3) (fcifc2 + fe2fc3 + fc3fcl) 
Ct (fci + fc2 + fc3) 
(fci + fc2 + fc3)'^ + 15fcifc2fc3 
4 (fci + fc2 + fc"3)' 
(fci -I- fc2 -I- fc3) (fclfc2 + fc2fc3 + fc3fcl) 



-^24 



4 (fci + fc2 + fc3) 



(4.7) 



We see that the magnitude of the bispectrum depends 
on C+ which in turn depends on the new energy scale 
M*. Since we've found in jTF that the general relativis- 
tic value of the power spectrum for tensor perturbations 
is obtainable in the HMT framework when H ^ Af* , we 
can have the same conclusion as that for scalar pertur- 
bations presented in [1^: a large bispectrum is possible, 
provided that Af* is not too much lower than the Planck 
scale. 



V. SHAPE OF THE BISPECTRUM 

We are now ready to plot the shapes of the bispec- 
trum. For si — S2 = S3 ~ 1 and si = S2 = —S3 ~ 1, we 
plot the shapes contributed by various terms separately 
in Fig.'s[T]and[2J These two figures represent all possi- 
ble configurations when we do not have parity violating 
terms. 

When spins of all 3 tensor fields are the same (si = 
S2 = S3 — 1), the signal generated by the relativistic 
term R peaks at the squeezed limit (fcs/fci — ?► 0), while 
the term KjElRf favors the equilateral shape (fc2/fci — 
k^/ki ~ 1). It's interesting to note that the other two 
terms, K^^Rij and (V*i?^'^) {\7iRjk), also generate larger 
signal in the squeezed limit, though they are of higher- 
order derivatives. This is indeed expected, if one realizes 
that the fc-dependence of them are similar to that of the 
GR term [cf. Eqs.( [3l2pl5| ]. 

When the spins are mixed (si = S2 = —S3 = 1), the 
terms R, R^^ Rij and {y^R^'^) iy iRjk) generate shapes 
similar to the previous case. A particular interesting re- 
sult is that the signal generated by R^.R'^R'^ no longer 
favors the equilateral shape but peaks in between the 
equilateral and squeezed limits. This can be understood 
as that for the mixed spin the product of the polariza- 
tion tensors gives a strong favor of the squeezed shape. 
To illustrate this, we plot the fc-dependence of them in 
both cases (-1- -|- -t-) and (-1- H — ) in Figs. Islandlij where 
we've defined 
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ConfLguration4 



Ci - -C2 - -C3 






2;,2 
3 



^1^2 



c. 



Ci - -C2 
Ci - -C3 



(123)+ cyclic, 
(123)+ cyclic, 



(123)+ cyclic. 



(123) + cyclic. 



(5.1) 






VI. CONSTRAINTS FROM PLANCK 
OBSERVATIONS 

Finally, we would like to comment on the recently- 
released results of Planck on primordial non-Gaussianity 
[21] . As noted at the end of Section IV, the magnitude of 
the bispectrum is dependent on the mass scale M* . Here 
we use the Planck result to obtain a constraint on Af* . 

Very roughly speaking, the non-linearity parameter for 
bispectrum can be estimated as 



(hhh) 
{hh){hh) 



{hhh) 



(5.1) 



where A^ is the scale-invariant power spectrum of the 
tensor perturbations. It is natural to assume that the 




FIG. 1: Shapes of (fcifc2fc3)^^G+++ (K) contributed by vari- 
ous terms. All are normalized to unity for equilateral limit. 
X = k2/ki,y = fcs/fci. 





(a) Configuration 1 




(b) Configuration 2 




(c)Configuration 3 




(d) Configuration 4 

FIG. 3: Shapes of the difTerent configurations of the polariza- 
tion tensors for si — S2 — S3 = 1. All are normalized to unity 
for equilateral limit, x = ^2/^1,3/ ~ k-i/k\. 



FIG. 2: Shapes of (feifefcs) ^G+h (K) contributed by var- 
ious terms. All are normalized to unity for equilateral limit. 
X = k2/ki,y = ks/ki. 




(a) Configuration 1 




(b)Configuration 2 




(c) Configuration 3 




(d) Configuration 4 

FIG. 4: Shapes of the different configurations of the polariza- 
tion tensors for Si = S2 — —S3 = 1. All are normalized to 
unity for equilateral limit, x = k2/k\,y = ka/ki. 



bispectrum of the tensor perturbations has a lower mag- 
nitude than does the scalar perturbations, 



Hence 



(hhh) < {nnn). 



{hhh) ^ _2{Tzmz) 



(5.2) 



-2 r-R 



r7Z 



(5.3) 



where r is the tensor-to-scalar ratio and A^ is the scale- 
invariant power spectrum of the scalar perturbations 
which takes the value (4.9 x 10^^) . 

The equilateral signal in our model takes the approxi- 
mate value 



(hhh) 



HC4 



96, 



(5.4) 



where we've taken il/pi = 1. The non-linearity parameter 
for the equilateral shape is constrained by the Planck 
result /^L ^ -42 ± 75, or |/^J ^ 100. This, together 
with Eqs.(5.3) and (5.4), requires 



HC. 



1.961 < 100 (4.9 X 10" 



-5^4 < 



(10-)' 



(5.5) 



The normalization condition for the mode function vj^ 

(5.6) 



in the relativistic region in (4.3 1 reads 
Vkv'C - vlv'k = ih, 
which requires 

^- . 1. 



(5.7) 



In the mean time, good IR behavior (so that the theory 
flows to GR in this limit) requires Cj, — I, hence C+ ~ 1. 
This implies 



Ct 



1 



a 



Apply this to condition (5.5), we have 
H'\96\ < (10-^)^ 



(5.8) 



(5.9) 



Since g^ and gg are both coupling constants for the cubic 
curvature terms, it's natural to assume that 



156 1 



This gives us the final constraint 



M* 



M, 



3/2 

EffL 



AL 



pi 



M, 



pi 



H 

Ml 



5,10" 



where we've written Mpi explicitly and e^i 
was introduced in 1181. 



(5.10) 
(5.11) 
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In the nonprojectable case without the extra U(l) 
symmetry, the frame effects imposed the most stringent 
constraint on the upper bound, M* < 10^^ GeV [S], 
while in the case with the local U(l) symmetry, such 
effects have not been worked out yet, either with the 
projectable condition or without it. Assuming that this 
is also true in the present case, the above condition im- 
plies that H/KU < 10"^. For M* ~ M^i, we find that 
H/M^ < 10-3. 

On the other hand, in jT8j it was found that when 
ehl = {H/M^^Y "^ 1) our model can reproduce the power 
spectrum given by the inflationary models i n GR . Tak- 
ing chl ^ 0(10^^), we find that constraint (5.111 gives 
M^/Mpy < 10-^ 



VII. CONCLUSIONS AND REMARKS 

In this paper, we study 3-point correlation function of 
primordial gravitational waves generated during the de 
Sitter expansion of the universe in the framework of the 
general covariant Hofava-Lifshitz gravity with the pro- 
jectability condition and an arbitrary coupling constant 
A. We find that the interaction Hamiltonian, at leading 
order, receives contribution from four terms built of the 
3-dimensional curvature: 



R, 



R^JR'', 



RjRiR^ 



(V'i?J''=) (V,i?,fc) . 



The Ricci scalar R yields the same fc-dependence as that 
in general relativity, i.e. its signal peaks at the squeezed 
limit regardless of the spins of the tensor fields, but with 
different magnitude due to coupling with the U{1) field A 



and a UV history when the dispersion relation is signifi- 
cantly different from the relativistic form. The two terms 
RijR^^ and (V'i?-''^) (ViRjk) generate shapes similarly 

to the R term. The term R^R-j^R^ favors the equilateral 
shape when spins of the three tensor fields are the same 
but peaks in between the equilateral and squeezed limits 
when spins are mixed. We find that this is due to the 
effect of the polarization tensors: when spins are mixed, 
the product of the three polarization tensors strongly fa- 
vors the squeezed shape. Thus, an absence of the equilat- 
eral shape in the bispectrum of the tensor perturbations 
cannot rule out gravity theories of higher order deriva- 
tives, at least for the R^ type of theories. 

It should be noted that in performing the integration 
over the three mode functions to get the full expression 
of the bispectrum [cf. Eq.(4.5)], we integrated over only 
the region when k^ dominated the dispersion relation [cf. 
Eq.(4.3)]. We gave a qualitative argument for the condi- 



tion under which the UV history can be (partly) ignored, 
in leading order analysis. A quantitative study of the er- 
rors introduced with such ignorance will certainly deserve 
further analysis. 

Finally, we've also obtained a constraint on M* and H 
[cf. Eq.(5.11)] using the Planck results, released recently 
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